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I. INTRODUCTION. 

Charmonium light cone wave functions (LCWF) are universal nonperturbative objects that describe the production 
of charmonium mesons in hard exclusive processes within light cone formalism 1]. Usually to study hard exclusive 
processes with quarkonium production one uses NRQCD 2]. So, light cone formalism can be considered as alternative 
to NRQCD. 

There are two very important advantages of light cone formalism in comparison to NRQCD. The first one is 
connected with the following fact: light cone formalism can be applied to study the production of any meson. For 
instance, it is possible to study the production of light mesons, such as n mesons, or the production of heavy mesons, 
such as charmonium mesons, if LCWFs of these mesons are known. From the NRQCD perspective, this implies that 
light cone formalism resums infinite series of relativistic corrections to the amplitude, what can be very important 
U S B S S The second advantage is that light cone formalism easily resums leading logarithmic radiative 
i— — corrections to the amplitude ~ a s log(Q) with the help of LCWFs. This is very important advantage since leading 

' logarithmic corrections at high energies can be even more important than relativistic corrections to the amplitude. 
Q-f From this one can conclude that LCWFs are the key ingredient of light cone formalism. Moreover, the universality 
of LCWFs and the variety of the processes where these functions can be used make the study of charmonium LCWFs 
to be a very important task. However, despite the fact that charmonium LCWFs are very important in understanding 
tIh ■ hard exclusive processes with charmonium production there is a very limited knowledge of the properties of these 
~~ '■ functions. There are only few papers where this functions were studied [H, ITol [Til. [l2| . 

\ In this paper the procedure developed in papers [H, [T3| for the study of charmonium LCWFs will be applied to 
1 the study of leading twist LCWFs of and rj' c mesons. This paper is organized as follows. In the next section all 
[ definitions needed in the calculation will be given. In Section III the moments of LCWFs will be calculated in the 
OO . framework of Buchmuller-Tye and Cornell potential models. Section IV is devoted to the calculation of the moments 
within NRQCD. QCD sum rules will be applied to the calculation of the moments in Section V. Using the results 
obtained in Sections III-V the models of LCWFs will be built in Section VI. In the last section the results of this 
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II. DEFINITIONS. 



There is one leading twist light cone wave function (LCWF) of r( c meson (f>ri(£,, \i) and there are two leading twist 
LCWFs of tp' meson /x), 0t (£,/•*)• The function <^j,(£,/z) is twist two LCWF of longitudinally polarized ip' 

meson. The function <^x(£, fi) is twist two LCWF of transversely polarized t/j' meson. These LCWFs can be defined 
as follows [l| 

(0\Q(z) lal5 [z,-z]Q(-z)W c (p)), = ifvPcJ dteWe^fon), 
(0\Q(z) la [z^z}Q(-z)\^(e x=0 ,p))^ = f L p a J ^ d£ e^^(£, fi), 
(0\Q(z)<j a p[z,-z}Q(-z)W(e x=±1 ,p))^ = Mn){e*pp ~ e?P«) J dte^fafon), (1) 

where the following designations are used: Xi,X2 are the momentum fractions of the whole meson carried by quark 
and antiquark correspondingly, £ = x\ — p is the momentum of corresponding meson, /i is an energy scale. The 
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FIG. 1: The functions $(£) for 15 and 25 states. 



factor [2,-2], makes the matrix elements to be gauge invariant and the dependence of the LCWFs <Po,l,t(x, A 4 ) on 
scale /x can be found in p], [ID, [l4| • 

It should be noted here that there is an important distinction between LCWF of rj' c and ip' mesons. Let us, for 
instance, consider LCWF of longitudinally polarized ip' . Obviously, this function can be written as follows 

^(£,m) = <*!(£,aO + (2) 

where </>£(£,//) and (£, /z) are S- and D-wave contributions to LCWF of ip' meson. In the case of r]' c meson only 
S'-wave contributes to the LCWF of this meson. It is not difficult to estimate the contribution of -D-wave to LCWF 
^(CjM)- Evidently, D-wave admixture in the LCWF (/>l(£, A*) is proportional to the factor ~ tan(6*)/f//f , where 
fhi fE are &~ an< i -D-wave contributions to the constant Sl, 9 is a mixing angle of S and D waves in ip' meson. Within 
potential models this factor can be written as 

tan(fl)/f I Si = tan(fl) J 3 R '^°} (3) 

where M c is a quark mass in the framework of potential model, Rg(r), Ru{r) are radial wave function of D and S waves. 
Numerical values of parameters 8, Rp(0), Rg(0), M c needed for the estimation of D-wave contribution to LCWF of 
iP' meson will be taken from paper 0: R s (0) = 0.734 GeV 3 / 2 , 5i?^,(0)/V8Af c 2 = 0.095 GeV 3/2 , 9 - 12°. Thus one 
gets rather large suppression of D-wave admixture tan(6')/f 5 //£ ~ 0.03. On account of the considerable suppression 
D-wave admixture one can disregard its contribution to the LCWFs of ip' meson. Below this approximation will be 
used. 

Commonly, r( c and ip 1 mesons are considered as a nonrelativistic bound states of quark-antiquark pair. At leading 
order approximation in relative velocity of quark-antiquark pair rj' c and ip' mesons cannot be distinguished. So within 
this approximation r]' c and ip' mesons have identical LCWFs at scale ~ M c 

M) = <M£,Ai) = <M£>M) = 4>(t,»)- (4) 

One can expect that in the case of 2S* mesons corrections to this approximation can be large. However, the accuracy 
obtained in this paper does not allow one to distinguish LCWFs <prj t L,T( x i /•*)• F° r this reason approximation ([4]) will 
be used in this paper. 

The main goal of this paper is to calculate the LCWFs <P v ,l,t(£,, m) °f "0' an d rf c mesons. These LCWFs will be 
parameterized by their moments (£^,l,t)p at some scale. It is worth noting that, the LCWFs (JTJ) are £-even, so only 
even moments should be calculated. 



III. THE MOMENTS IN THE FRAMEWORK OF POTENTIAL MODELS. 

In papers [TH, [3| it was shown that the moments of LCWFs of i] c and J/ip mesons can be calculated in the 
framework of potential models. In comparison with QCD sum rules, such calculation cannot be considered as an 
accurate one. However, potential models give rather good estimation of the values of the moments. 

To calculate the moment of LCWF one can apply Brodsky-Huang-Lepage (BHL) [l6| procedure that can be written 

as 

~ as (O<i ) (^M) = (i-e 2 )$(^/i), 
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TABLE I: The moments of LCWF obtained within different approaches. In the second and third columns the moments 
calculated in the framework of Buchmuller-Tye and Cornell potential models are presented. In the fourth column NRQCD 
predictions for the moments are presented. In last column the results obtained within QCD sum rules are shown. 



where ip(k 2 ) is the solution of Schrodinger equation in momentum space, M c is a quark mass within potential model. 

In this paper the function ~0( k 2 ) will be calculated in the framework of the potential models with Buchmuller-Tye 
[r^ and Cornell potentials [la ]. The parameters of Buchmuller-Tye potential model will be taken from paper p7| . 
For Cornell potential V(r) = —k/r + r/a 2 the calculation will be carried out with the following set of parameters: 
k = 0.358, a = 2.381 GeV" 1 , M c = 1.147 GeV [9]. The scale n is taken equal to 1.5 GeV. 

The results of our calculation are presented in Table I. In second and third columns the moments calculated in 
the framework of Buchmuller-Tye and Cornell models are presented. It is seen that there is good agreement between 
these two models. In papers [13, 0] it was shown that potential models cannot be applied for higher moments. Due 
to this fact the calculations have been restricted by few first moments. 

Now the following point deserves consideration. It was shown in equations ([5|) that LCWFs can be represented as a 
product of the asymptotic function as (£) and the function <!>(£,//). The function <!>(£, /i) contains information about 
the internal motion of quark antiquark pair in meson. Let us compare these functions $(£) for IS" and 2S states. 
These functions calculated within Buchmuller-Tye potential model and normalized as J = 1 are shown in Fig. 

1. It is seen that the function <&(£) of 15* state has rather simple shape with one extremum at £ = 0. It is not difficult 
to guess that this extremum appears since the function $(£) is £— even. For the same reason the function $(£) of 2S 
state has similar extremum at £ = 0. However, the function of 2S state has two additional extremums located 
symmetrically relative to £— axis. 

It is not difficult to understand why these additional extremums appear. To do this let us differentiate equation 
© over £ 



(M 2 + fi ) ip \ c ^ - M 2 C V 



i - f 2 / V 1 - £ 



(6) 



Equation ([6]) can be simplified if one recalls that the scale /i is much greater than characteristic momentum of relative 
motion of quark-antiquark pair inside the meson. This means that the function ip((£ 2 M 2 + /i 2 )/(l — £ 2 )) in the first 
term is much less than ?/>((£ 2 M 2 )/(l — £ 2 )) in the second term of equation ([6]) for not too large £. So, the first term 
gives small correction to the second and can be omitted to the first approximation. Then equation ([6]) can be written 
as 



It is well known that equal time wave function i/^k 2 ) of 2S state has one zero at some point k 2 . So it is clear that 
the function <&'(£, fj) changes sign at the points £ = ±ykjj/ (k 2 ± M 2 ) what corresponds to the two extremums of the 
function <!>(£, n). Moreover, the function $'(£, n) changes sign at £ = what corresponds to the extremum at £ = 0. 
Obviously, if one regards the first term in equation © this will just shift the position of extremums. 

Applying the same arguments it is not difficult to prove the following statement: leading twist LCWF of nS 
state has 2n+l extremums. It should be noted here that our arguments are based on the relation between LCWF 
of leading twist and equal time wave function ((5]). In papers 0, [T3| the other relations were proposed . Nevertheless, 
the statement written above remains true since it is valid for all relations of the type (f51) where the function ip(t) can 
be represented as a product of equal time wave function and some function x(C) ~ 1 + 0(v 2 ) for £ ~ v. 
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FIG. 2: Sum rules for with different values of parameter fig. a (££}=0; fig. b <£Z}=0.22; fig. c <££)=0.18. 



IV. THE MOMENTS IN THE FRAMEWORK OF NRQCD. 



To calculate the moments of LCWFs at leading order approximation in relative velocity one can use the following 
formula (l3j : 

71 + 1 

where the constant 7 can be related to the matrix element of NRQCD operator j 2 — (v 2 ). The value of (v 2 ) can be 
calculated using the approach proposed in 

(v 2 ) = 0.65 ± 0.42. (9) 

There are different sources of error to result ([9]). However, the main source of error is relativistic corrections to 
formula ©.In © the size of these corrections was estimated as ~ {(v 2 )) 2 . It is interesting to note that the value 
© obtained at leading order approximation in relative velocity is very close to that obtained at next to leading order 
approximation (v 2 ) = 0.67 [201 ] . 

Using ([8]), ([9]) one can easily calculate the values of the moments. The results of this calculation are presented in 
the fourth column of Table I. The central values of the moments were calculated according to formulas ([8]). The errors 
of the calculation of the moment (£ 2fc ) were estimated as ~ k(v 2 ) x (£ 2fe ). 

It is seen from Table I that within the error NRQCD prediction for the second moment is in agreement with potential 
model estimation, but the central values are rather far from each other. For higher moments the difference between 
central values obtained within these approaches becomes more dramatic and the errors of the calculation within 
NRQCD are very large. From this one can draw a conclusion: although NRQCD can be applied to the calculation 
of the second moment of 2S state mesons, the predictions obtained within this approach for higher moments become 
unreliable due to large relativistic corrections. 

It should be noted here that formula © is very simple. So it is not difficult to guess that this dependence can be 
reproduced by the following function 

0(0 = ^0(7-1*1), (10) 

where 9(x) is the Heaviside step function. Function (JTUJ) can be considered as the NRQCD LCWF obtained at leading 
order in relative velocity. If the velocity of quark-anitquark pair is infinitely small (7 — > 0) than <fr(£) tends to as 
it should be. 

Function (|10[) is very simple and it does not reproduce peculiarities of mesons. For instance, the only distinction of 
LCWFs of IS and 25 states is different constants 7 of these mesons. However, from consideration of previous section 
it is known that the forms of these LCWFs are rather different. Actually, this is not surprising if one recalls that 
within NRQCD all mesons with the same quantum numbers are described identically by one set of constants and the 
peculiarities of each meson are contained in the values of these constants. At leading order approximation in relative 
velocity there is only one constant (v 2 ). So it is not possible to reproduce peculiarities of each meson by the only 
constant. Probably, if one regards relativistic corrections and QCD radiative corrections to the expressions ((SJ some 
properties will be restored. 
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V. THE MOMENTS IN THE FRAMEWORK OF QCD SUM RULES. 
A. The moments of /i). 

In this section QCD sum rules [HI, [22j will be applied to the calculation of the moments [H, [23| of LCWFs </>l(£, /i). 
To do this let us consider two-point correlator: 

U L (z,q,n) = i f dSe^<0|TJ (aOJ„(0)|0) = (zq) n+2 U L (q 2 , n), (11) 



J (x) = Q{x)zQ{x) 7 J„(0) = Q{0)z{iz p D p ) n Q{0), z 2 = 0. 
Sum rules for this correlator can be written as follows: 

(Zl) j/* (££).//* (Jx)|'(£lV i f so imn pcrt (s,n) ( m ) 2 2 

(Mj^ + Q 2 ) m + l + (M 2 ,+Q^ = n L, dS (s + Q^ + > ») = • »)• ^ 

where The expressions for the functions Imn pcrt (s,n) and n| 1 ™] rt (Q 2 , n) can be found in paper [H]. (/l),//* and 
(/i)^' are leptonic constants of J/^f and -0' meson, and (C£)v are the n-th moment of J/ 1 ^ and mesons' 

LCWFs. To remain the designations introduced earlier, below and (££) will be used instead of ( fr. )■,/, > and (££)^/>'- 
Numerical analysis of QCD sum rules (|12j) will be done similar to the numerical analysis in paper [1 31 ] - To weaken 
the role of unknown radiative corrections instead of sum rules (fT2j) the ratio of sum rules with different n will be 
considered: 

(&)j/* + r (&)a(m) _ n L (Q 2 ,0) 



l + ra(m) II L (Q 2 ,n)' 



where 



r = /!7(.fc) 2 //*, 



ttH = f-^±gr +1 , a w^c+o * o.53 d4) 

To calculate the moments of LCWF (/>l(£,,I^) let us rewrite sum rules (fT3"|) as 

(S>J/ * = S(QV) (1 + r ° (W)) " <S> " fl(m) ' (15) 

First sum rules (|15p for n = 2 will be considered. To the first approximation let us disregard the contribution of 
meson in the right hand side of equation (fTS)) . as it was done in paper [lj] and take the value of the threshold so 
equal to the threshold of D-mesons production y/s ~ 3.7 GeV. The left hand side of equation (fT5|) does not depend 
on to. The right hand side of (fl5|) is a function of to. This function is plotted in Fig. 2a. It is seen that for too 
small values of to (to < 10) right hand side of equation (|15[) varies rather rapidly. This happens since there are large 
contributions from higher resonances disregarded in model of physical spectral density what invalidates sum rules 
(p~2|) . (fTS"|) . Although for to 3> m\ these contributions are strongly suppressed, it is not possible to apply sum rules 
for too large to (to > 12) since the contribution arising from higher dimensional vacuum condensates rapidly grows 
with m,(see Fig. 2a) what also invalidates sum rules. It is seen from Fig. 2a that in the region [10, 12] left hand side 
of equation (fT5|) m varies very slowly. This is the region of applicability of sum rules (fT2"]) . (fT5|) where the resonance 
and the higher dimensional vacuum condensates contributions are not too large. Within the region of applicability 
the approximation of physical spectral density and the approximation of the contribution of vacuum condensates are 
valid and one can determine the value of the constant (C£) j/*- Thus one gets 

= 0.07. (16) 

This value coincides with that found in paper jl4j |. 

As it was noted above due to the contribution of higher resonances sum rules (|15p is spoiled in low to region. 
Evidently, the inclusion one resonance succeeding J/^-meson will improve sum rules (115[) in the region of low to. The 
parameter (£ 2 ) can be chosen so that to attain best fit of right hand side of equation (| 1 5|) to the constant (££),//* . The 
calculation shows that the best fit can be obtained if (£|) = 0.22. Right hand side of sum rules (fT5|) at (£ 2 ) = 0.22 
as a function of to is shown in Fig. 2b. 
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FIG. 3: Allowed region for parameters (/3, a) ( model (|21[1 ) is painted black. 



From Fig. 2b it is seen that if ip' meson with = 0.22 is included into the sum rules, the agreement between right 
and left hand sides of equation (fT"5|) becomes much better. From Fig. 2b one also sees that in the region m e [0, 4] right 
hand side of sum rules (fTS"|) is rising function of m. This seems rather strange since if one includes charmonium meson 
succeeding ip' meson to sum rules, right hand side of equation (| 1 5 j) will become decreasing function of m. Perhaps, this 
strange behavior originates from the following fact. In the region of too low m there are large contributions coming 
from higher resonances not included into physical spectral density. So, if one tries to regard these contributions by 
the only resonance - ip' meson, this will lead to an overestimation of the value of (££)■ This problem can be partially 
removed if, in addition to the requirement to achieve the best fit of both sides of sum rules, the following requirement 
will be imposed: right hand side of equation (TT5)) must be decreasing function of m. Thus one gets = 0.18. The 
right hand side of sum rules (fl~5|) as a function of m with = 0.18 is plotted in Fig. 2c. 

There are many sources of uncertainty of the calculation fulfilled above. The first one appears due to the uncertainty 
in sum rules parameters m c and (a s G 2 /7r) (l4| . The calculation shows that the uncertainties due to the variation of 
m c and {a s G 2 /it) are not very important (not greater than 10%). For this reason this source of uncertainty will not 
be considered in the calculation. Probably, the unknown contribution of QCD radiative corrections to the spectral 
density is much more important, but it is difficult to estimate its contribution. Another very important source of 
uncertainty results from the unknown value of the threshold parameter \/s Q - This parameter determines the energy 
from which continium contribution to sum rules appears. It is difficult to calculate the value of so, one can only claim 
that it is not very far from the threshold of D-mesons production ^/s ~ 3.7 GeV. In the calculation carried out 
in this paper it will be assumed that ^/s Q belongs to the interval 3.7 ± 0.5 GeV. The interval chosen in such a way 
is rather broad and it contains all intervals common for QCD sum rules analysis. It should be noted here that the 
error due to the variation of sq within this interval is rather large and below it will be considered as the error of the 
calculation. 

Applying the method discussed above for higher moments one gets the results: 



The central values of the moments have been calculated at y/s — 3.7 GeV. The errors of the calculation appears 
due to the variation of the threshold parameter y/s within the interval 3.7 ± 0.5 GeV. Physically this variation can 
be considered as a simulation of the contributions of higher charmonium states and continuum to the moments of 
LCWF </>l(£, m)- From this perspective the error of the calculation is rather large since the contributions from ip' 
meson, higher resonances and continuum are not well separated in sum rules (|11[) . All these contributions appear 
approximately at \fs — 3.7 GeV. So one can conclude that, although this source of uncertainty can be diminished, 
it will remain to be the main source of uncertainty of the calculation. From results (|17|) one sees that the error of 
the calculation rises as number of the moment increases. Evidently, this happens since the larger the number of the 
moment the larger the sensitivity of this moment to higher charmonium states and continuum. 

Results of the calculation (jTTJ) are presented in the fifth column of Table I. It is seen from this table that, although 
the accuracy of the results obtained within sum rules is better than NRQCD predictions for the moments, the error 
of the calculation is still rather large. It should be noted also that QCD sum rules predictions for the moments are in 
better agreement with potential models than with NRQCD results. The central values of NRQCD predictions seems 
to be overestimated. 
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FIG. 4: The LCWF JUJ at scales /i = 1.2 GeV,Mi = 10 GeV,At2 = 100 GeV,^ 3 = oo. 



B. The moments of $t(£,IJ<) and </>,, c (£, p). 

It is not difficult to derive sum rules for 0t (£,/•*) and ^(Ca 4 )- F° r instance, to calculate the moments of ^(CjA 4 ) 
one should consider two-point correlator: 

n,(«,9,n) = i J d*xe* x {0\TJ Q (x)J n (0)\0) = (zq) n+2 U v (q 2 ,n), 

J (x)=Q(x) l5 zQ(x), J n (0)=Q{0) l5 z(iz p D p ) n Q(0), z 2 =0. (18) 
Sum rules for this correlator can be written as 

fl^h | , 1 r . Imn pert ( S ,n) M 

(M2 + Q 2 )™+1 + (M| ol + Q 2 )"^ 1 + (M2 + Q2)m+1 ^ 7^ ( s + Q2 )m +1 f npcrtW ' 

where (£, n ) Vc , (£ n )xd and (£™/V are moments of leading twist LCWF of r; c , Xcij ?7e mesons, the constants / 7)c , / Xol , ,/y 
are defined as 

(0|Q(0) 7q75 Q(0)|M(p)) - if M Pa, M = r, c , X ci,v'c- (20) 

One sees that in addition to r]' c meson there is contribution of Xd meson. Since Xd meson is P wave meson, its 
contribution is a little bit suppressed. Nevertheless, sum rules (| 19[) has one additional unknown parameter (£ n ) Xcl 
and this leads to worsening of sum rules predictions in comparison to case considered above. Similar situation takes 
place for 4>t{£,t^), where there is contribution of h c charmonium meson. 

From this one can conclude that, unfortunately, QCD sum rules cannot distinguish LCWF 0l(£,/x), </>t(£, A*) and 
4>rj c (CiA*) and it is not possible to calculate the moments of 0t(£,A*) and ^ c (^, /x) with the accuracy better than the 
accuracy of the moments (£f,)j/*. This makes the calculation of the moments of </>t(£, A*) ^rj c (£,A*) within QCD sum 
rules rather pointless. Below hypothesis ^ with moments (fT?)) will be used. 



VI. THE MODEL FOR THE FUNCTIONS ^, l ,t(x,h). 



Unfortunately, the methods applied in this paper to the calculation of the moments do not allow one to distinguish 
LCWFs 4>n,L,T{x, /i)- For this reason, below these functions are assumed to be equal to some function <p(x, fj,) at scale 
fi ~ m c . This section is devoted to the construction of the model for this function based on the results obtained 
within QCD sum rules. Results (fT7| are defined at scale [i ~ m c [13]. In the calculations it will be assumed that 
these results are defined at scale fio = 1.2 GeV^ m c . 

In papers [13, 03] it was proposed one parametric model of LCWFs of r\ c and J/^f mesons at scale (j.q = 1-2 GeV. 
To reproduce the results obtained in this paper this function can be modified by additional factor (a + £ 2 ) 

0(6/* = A*o) = c(a,/3)(l - £ 2 )(a + £ 2 )exp(- 7^2) = <<*,P)Q- - = Mo), (21) 

Potential model calculation of <E'(£,A i ~ m c) tells us that this function is positive and it has three extremums. Below 
it will be assumed that these properties remain true for real function <&(£, /1 = a*o)- To meet the first requirement 
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<H£, /xo = 1.2 GeV) 


<?!>(£, Mi = 10 GeV) 


4>{£,H2 = 100 GeV) 


<K£,M3 = oo) 


n = 2 


0.18 


0.19 


0.19 


0.20 


n = 4 


0.051 


0.068 


0.074 


0.086 


n = 6 


0.018 


0.032 


0.037 


0.048 



TABLE II: The moments of LCWF proposed in this paper at scales fi = 1.2 GeV,/Xi = 10 GeV,/x 2 = 100 GeV, tt 3 = oo 
are presented in second, third, fourth and fifth columns. 

one can suppose that a > 0. It is not difficult to show that the function // ~ m c ) has extremums at £ = 0,£ 2 = 
(2 + a — ^/(2 + a 2 ) — 4(1 — a(3)) /2. Two additional extremums of the function /u ~ m c ) are beyond the physical 
region £ G [—1, 1]. So, to meet the second requirement - the function /x = fJ-o) must have three extremums - one 
should impose the condition a/3 < 1. 

Further let us find the region where the constant (3 can vary. This can be done in the framework of Borel version of 
QCD sum rules [24( where this constant can be expressed through the Borel parameter M as follows (3 — Am^/M 2 . 
The value of Borel parameter cannot be too small (M > 1 GeV), otherwise the vacuum condensates contributions 
become too large. At the same time Borel parameter cannot be too large (M < 3 GeV) otherwise the contributions 
of higher resonances become too large. Thus one gets the assessment of the interval where the constant f3 can vary 
(3 € (0.69, 6.25). Now it causes no difficulties to find allowed region of the constants a, (3. This region is painted black 
in Fig. 3. 

The central values of the second and the forth moment can be obtained within model (|2~lj) if the values of the 
constants {a, (3) are equal to (0.027,2.49). If one fixes the value of the constant a = 0.027 than, to attain the 
agreement of the model ([2~Tj) with the results (fl7| for the second moment, the constant (3 can vary within the interval 
(3 e (1.4, 5.7). Similarly if the constant (3 is fixed at 2.49 than the constant a can vary within the interval a £ (0, 0.35). 

Now let us consider model (f2~Tj) with the central values a — 0.027, (3 = 2.49. As it was noted above model (f2"Tj) with 
these values of the constants a, (3 is defined at scale fi = (j,q. It is not difficult to calculate this function at any scale 
fi > /Lto using conformal expansion ; lj. This calculation will be done only for the function </fi(x,/i). The function 
4>l{x,h) at scales /xo = 1.2 GeV, \i\ — 10 GeV, ^2 = 100 GeV, /X3 = 00 are shown in Fig. 4. The moments of this 
LCWF at scales (j,q — 1.2 GeV, /ii = 10 GeV, [i% = 100 GeV, /X3 = 00 are presented in second, third, fourth and fifth 
columns of Table II. 

In papers [HI, 03] it was shown that due to evolution LCWFs of IS* state have some interesting properties: the 
violation of nonrelativistic QCD velocity scaling rules, appearance of relativistic tail and improvement of the accuracy 
of the model. LCWFs of 25* states have similar properties and in this paper these properties will not considered. 

Now let us consider two different models jH}: Model I (a = 0,/3 = 2.5) and Model II (a = 0.2,0 = 2.5). LCWF 
<K£jM = Mo) °f these models are shown in Fig. 5a. LCWF of Model I has the following moments (£ 2 ) = 0.21, (£ 4 ) = 

0. 061, Model II has the moments (£ 2 ) = 0.12, (£ 4 ) = 0.031. It is seen that Model I is considerably wider than Model 
II. In addition, Models I and II are physically different. Really, suppose the meson with momentum p has LCWF 
of Model I. It is seen from Fig. 5a that this LCWF has rather sharp extremums at |£| ~ 0.5. This means that 
within this model it is not possible to produce 25* state charmonium meson from quark-antiquark pair with small 
relative momentum. Contrary to Model I, within Model II it is possible for quark-antiquar pair to have small relative 
momentum. Unfortunately, the uncertainties of results (|17l) are rather large. So, both models are allowed. One can 
only assert that the model of LCWF with central values of parameters a — 0.027, (3 = 2.49 is very similar to Model 

1. In addition, the forms of LCWF obtained within potential models (see Fig. 1) are similar to Model I. It should 
be noted here that at leading order approximation of NRQCD quark-antiquark pair has zero relative momentum. So 
this approximation is in contradiction with Model I. 

The effect considered above takes place at scale /i = /xq. To understand what happens at larger scales one should 
evolve Models I and II from scale /xo to larger scales. LCWFs of Models I and II at scale /x = 10 GeV are shown in 
Fig. 5b. It is seen from this plot that the effect is not so dramatic as it is at scale fj,Q. This result is in agreement with 
the property of LCWFs discussed above: the larger the scale the less difference between different models of LCWF. 

VII. CONCLUSION 

In this paper the moments of leading twist light cone wave functions (LCWF) of 2S state charmonium mesons have 
been calculated within three approaches. In the first approach Buchmuller-Tye and Cornell potential models were 
applied to the calculation of the moments of LCWFs. In the second approach the moments of LCWFs were calculated 
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in the framework of NRQCD. In the third approach the method QCD sum rules was applied to the calculation of 
the moments. Although, the results of the calculation are in reasonable agreement with each other, the errors of the 
calculation are rather large. As the result, it is not possible to distinguish different LCWFs form each other. 

Similarly to the study of LCWFs of IS" state charmonium mesons [HI, [HI, the most accurate results were obtained 
within QCD sum rules. Using these results two parametric model of LCWFs of 2S states was proposed. This model 
can be used in the calculation of different hard exclusive processes with 2S charmonium mesons production. 
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